Universality and nonmonotonic finite-size effects above the upper critical dimension 
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We analyze universal and nonuniversal finite-size effects of lattice systems in a L d geometry above 
the upper critical dimension d — 4 within the 0(n) symmetric ip 4 lattice theory. On the basis of 
exact results for n — > oo and one-loop results for n = 1 we identify significant lattice effects that 
cannot be explained by the ip 4 continuum theory. Our analysis resolves longstanding discrepancies 
between earlier asymptotic theories and Monte Carlo (MC) data for the five-dimensional Ising model 
of small size. We predict a nonmonotonic L dependence of the scaled susceptibility \L^ d ^ 2 at T c 
with a weak maximum that has not yet been detected by MC data. 
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The concept of universality plays a fundamental role 
in statistical and elementary particle physics . ft im- 
plies that a unifying description of various physically dif- 
ferent lattice and continuum systems near criticality can 
be given within the (p 4 field theory with the Hamiltonian 



H= d d 



ro 2 



(1) 



The wide applicability of this theory is well established 
below the upper critical dimension d* = 4 [0J§]. Partic- 
ular accuracy has been achieved in testing the universal 
predictions of the ip 4 theory by means of numerical data 
for the universality class of the d = 3 Ising model not only 
for bulk properties but also for finite-size effects with pe- 
riodic boundary conditions (p.b.c.) 

Less well established, however, is the range of appli- 
cability of the ip 4 theory for confined systems above the 
upper critical dimension where the critical exponents are 
mean- field like [jl|j2| . Early disagreements between Monte 
Carlo (MC) data for the finite d = 5 Ising model Q and 
universal predictions based on H jj] have led to a long- 
standing debate |j . New discrepancies between accurate 
MC data [Hi and recent quantitative finite-size scaling 
predictions [11 based on the (p 4 lattice Hamiltonian 



H 



E ofo 



^) 2 (2) 



<ij> 



have raised the question to what extent the <p 4 theory is 
capable of describing finite-size effects of the Ising model 
for d > 4. In particular the recently discovered |9|,|l0| non- 
equivalence of H and H for finite systems is in striking 
contrast to the situation for d < 4. This non-equivalence 
may be relevant not only for higher-dimensional finite 
systems but also for three-dimensional physical systems 
for which mean-field theory provides a good description, 
such as systems with long but finite range interactions 
[0] '. polymer mixtures near their critical point of unmix- 
ing jl^] , and systems with a tricritical point |l3| . 

In this Letter we resolve the existing discrepancies for 
d > 4 on the basis of exact results for the O(n) symmetric 



tp theory in the limit n — > oo and of one-loop results for 
n = 1. Our analysis of both H and H with a smooth and 
a sharp cutoff is not restricted to large L and allows us to 
specify the range of validity of universal finite-size scal- 
ing for p.b.c. in a L d geometry. We find, for p.b.c, that 
H with a smooth cutoff belongs to the same universality 
class as H whereas H with a sharp cutoff exhibits differ- 
ent nonuniversal finite-size effects. This implies that the 
lowest-mode prediction [Q of universal ratios at T c for 
d > 4 is indeed valid asymptotically for both the lattice 
(p 4 theory and the continuum (p 4 theory with a smooth 
cutoff. We demonstrate, however, that the existing MC 
data for the d = 5 Ising model of small size [||-|| are 
outside the asymptotic scaling regime and cannot be ex- 
plained by H because of significant lattice effects. We 
also demonstrate that our one-loop results based on H 
are in quantitative agreement with the MC data || for 
4 < L < 22 and that the one-loop two-variable scaling 
results ||,|l4| are well applicable to L > 12, contrary to 
earlier conclusions J^Jlq]. We predict a weak maximum 
of the L-dependence of the scaled susceptibility xL~ d / 2 
at T c which has not yet been detected in the MC data M . 
Our analysis implies £o = 0.396 for the bulk correlation- 
length amplitude of the d = 5 Ising model, in disagree- 
ment with Co = 0.549 found in Ref. §. 

We start from H, Eq.(||), for the n-component vari- 
ables ipi on a finite sc lattice of volume L d with a nearest- 
neighbor coupling J > 0. The basic question is to 
what extent H is equivalent to the spin Hamiltonian 
H s = —K J2<ij> s i s j where the n-component spin vari- 
ables have a fixed length sf = n, in contrast to ipi whose 
components (pi a vary in the range — oo < tpi a < oo. For 
n = 1,H S is the Ising Hamiltonian with Sj = ± 1 and 
K>0. 

An exact equivalence between H and H s exists in the 
limit uq — > oo, ro — ► — oo at fixed uo/(Jro) for general L, 
n and d. Choosing uo/(Jro) such that K = —Jro/^iuon) 
we obtain by means of a saddle-point integration 



hm Y 

UQ-.QO 



K 
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Xs 



(3) 
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where x and Xs are the susceptibilities 

X = (nL*)- 1 Y,<ViVj>, (4) 
Xs = {nL d )- x <s lSj > . (5) 



h3 



The weights in Eqs. 



and fj) 



are e 



-H 



and 



-H, 



re- 



spectively. For n = 1, this exact equivalence is of limited 
relevance since all calculations within the (p 4 model are 
performed at finite uo. Hence, even in an exact theory, we 
have Xs ^ Jx/K at finite uq. Therefore, in a quantita- 
tive comparison of x with MC data for Xs , one must allow 
for a (T and L independent) overall amplitude A which 
is adjusted such that Xs — AJx/K. For finite u , the 
constant A accounts for an appropriate normalization of 
the variables ifi relative to the discrete variables Si = ±1. 
In an approximate theory, the value of A depends on the 
approximations made for x- This corresponds to an ad- 
justment merely of the nonuniversal bulk amplitude and 
not of the L dependence of x (f° r d = 3 see, e.g., Ref. 
Q). An adjustment of A was not taken into account in 
the analysis of Ref. || . 

Of particular interest is the case n — > oo since it pro- 
vides the opportunity of studying the exact Uq depen- 
dence including u$ — > oo. This reveals the structural 
similarity between x at finite uq and at uq — oo. This is 
most informative for d > 4 where the leading and sub- 
leading powers of L are independent of n and should 
apply also to the Ising universality class with n = 1. 

For n — > oo at fixed u^n the susceptibility x ~ 2 Jx f° r 
p.b.c. is determined implicitly by |10| 



X 



r /(2J) + J~ 2 u nL 



(0) 



with GkOT 1 ) = Or 1 + ^k)- 1 and J k = 2£i=i(i - 
cos kj) where ^ k runs over k vectors with components 
kj = 2Trrrij/L, nij = 0, ±1, ±2...,j = 1,2,..., din the 
range — 7T < kj < w. At T = T c we derive from Eq. (|^) 
the exact implicit equation for d > 4 

_ d AoM-x (4 - d)/2 Mr 1 ) 



x- 1 &i(x-\ L) 



(7) 



with Aq(uq) = (J 2 + uqtl J J k 2 ) (uon) 1 and 



Mr 1 ) = x {d - 6)/2 J [Jtir' + J^V 1 , (8) 

k 

A^rU) = J c k (r T 1 - L- d E c k (r T, w 



k^O 



where / k = (2ir)- d f d d k with \kj\ < n. We see that 
the structure of the L dependence of x f° r finite uq > 
is the same as for uq — > oo where Aq(mo) is reduced to 



Ao = Jk^k 2 - ^ i s reasonable to expect that also for 
n = 1 the calculation of x at finite Uq yields essentially 
the correct structure of Xs- 
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exact lattice theory, Eq.(7) 

scaling function, Ref. [10] 
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FIG. 1. Scaled susceptibilities for d — 5 at T c . Solid and 
dashed lines approach the lowest-mode lines for L —* oo. 

In Fig. la we show the exact result of %L~ 5 / 2 for 
n — > oo and d = 5 at T c by solving Eq. (]?]) numerically 
with A = J k J k 2 = 0.01935. We find that X L' 5/2 has a 
weak maximum at L = 9 which is not contained in the 
(large L) scaling form Xscai 

= L d ' 2 P{L 4 - d /X ) of Ref. 
p0{ (dashed curve). In Xscai the nonasymptotic Wegner 
correction oc was neglected and Ai was approximated 
only by the leading term Ai = 7i(x -1 L 2 ) L 2 ~ d with 



I m {x) 



t m ^[K b (t) d -K(t) d + l] 



where K b (t) = (ir/t) 1 / 2 and K(t) = E^-oo e M-ft)- 
Both x and Xscai show the predicted [§ slow 0{L {A - d V 2 ) 
approach to the large-L limit x^L^ d l 2 = X^ 2 corre- 
sponding to the lowest-mode approximation (horizontal 
line in Fig. la). Note that both x and Xscai approach xo 
from above. 

The small difference between x and Xscai in Fig. la for 
L > 15 arises from the negative Wegner correction term 
-X (4 " d)/2 /b(x _1 ) oc -L( 4 - d ) d / 4 / fc (0) in the numerator of 
Eq. (Q). The pronounced departure of x from Xscai for 
L < 10, however, is a lattice effect that is dominated by 
the subleading term —M\L~ d in 



A 1 (r 1 ,L) = h(x)L 2 



Mi{x)L- d + 0(L-"- 2 ), (11) 



-d-2\ 



2 



Mi (a;) 



dt 



[Kjt^K"^) - K^tY^K'W)] 

p {xt/4TT 2 ) 



(12) 



the order parameter distribution where $ = L d £j ipj. 
The analytic result reads for arbitrary L H 



with x = x L . Unlike the leading term 1\L , the 
lattice term —MiL~ d cannot be incorporated in the uni- 
versal finite-size scaling function P(y) which depends on 
y = (L/l ) 4 ~ d with l 4 ~ d = Xq. In summary, the leading 
L dependence of x is represented as 



X = L d ' 2 (u^)- 1 ' 2 tf 2 
Q = MY eff ) 2 / MY eff ), 

ye// =i <i/2 r e// (u e// ) _l/ 2j 

J^dss m exp(-^Ys 2 -s i ) 



1 



<12 



L (i-d)d/4 



1 - q 1 L( i - d )/ 2 + q 3 L- d / 2 



1/2 



j- dsex p(-ir s2 - s 4) 



(16) 
(17) 
(18) 

(19) 



(13) 



with the effective parameters 



where q x = X Q 1/2 Ii(x), q 2 = A a/4 f b (0), and q 3 = 

— 1/2 ^ 

A Mi (x) . The functions I\ (x) and Mi (x) have a weak 
x dependence with J x (0) = 0.107 and Mi(0) = 0.676 for 
d = 5. Eq. ( |l3| ) is shown in Fig. la as dot-dashed line 
which approximates the exact result, Eq. (0), with very 
good accuracy down to L = 3. 

Now we turn to the question to what extent H, Eq. ([!]), 
is equivalent to H. From our result of x, Eqs. (|)-(|), we 
obtain the corresponding result of x field = n^ 1 J d x < 
f(x)ip(0) > after replacing Jk by k 2 and setting 2 J = 1. 
A novel feature for d > 4 is the fact that Ai depends 
significantly on the cutoff procedure. We need to dis- 
tinguish two cases : (a) a sharp cutoff A which restricts 
the k vector to \kj\ < A, (b) a smooth cutoff A where 
— oo < kj < oo but where (x _1 + k 2 )~ m is replaced by 
the (Schwinger type) regularized form Q {x~ 1 + k 



-d/4 



r e ff = dot + 12u (Si - A ) + 144 u 2 Q M 2 S 2 , (20) 

(21) 

^ .. - „ - (22) 



= u - 36u S2, 



Sn 



(a i + 12u M 2 + J k )- 



Jl/n 



(L d ^ )- 1/2 ML d/2 dot u - 1/2 ) 



(23) 



/~2 ^ 



2^-m 
reg 



exp [— (x 1 + fc 2 )s]. The former case (a) 



implies §,0 Ai oc L 2 and x field ot L d 2 at T c which 
differs fundamentally from the lattice result X oc L d l 2 . 
In the latter case (b), however, Eqs. (|ll|) and (H) are 
replaced by 

AiOTU) = h(x)L 2 - d - M 1 {x 1 )L- d + 0(e- A2L2 ), (14) 
Mi(r 1 )=x[l-cxp(-r 1 A- 2 )], (15) 

with the same leading term I\L 2 ~ d . This implies that 
Xfieid with a smooth cutoff has the same asymptotic 
(large L) finite-size scaling behavior as Xscai- Adjust- 



field 



Jk(0 



reg 



to 



ment of the leading amplitude X 

the lattice counterpart Ao = f k J^ 2 fixes the cutoff as 
A = 0.185 and M 1 (Q) = A" 2 = 0.034 for d = 5 which is 
smaller than Mi(0) by a factor of 20. This difference be- 
tween Mi and Mi constitutes a significant lattice effect 
for small L that is exhibited in Fig. la, with Xfieid L~ 5 ^ 2 
represented by the dotted line. We conclude that H with 
a smooth cutoff yields the same (large L) finite-size scal- 
ing behavior as H (for cubic geometry and p.b.c.) but 
does not account for the strong L-dependence of xL~ d ^ 2 
for small L. We expect this conclusion to hold for general 
n. 

Now we consider H for the relevant case n = 1. We 
start from the one-loop result for x = 2Jx and for the 
ratio Q =< $ 2 > 2 / < $ 4 > of moments < $ m > for 



The r.h.s. of Eqs. @ - (||) depend only on the 
parameters Ho = ito/(4^P) and ao = ao/(2J) where 
a = (r - roc)/* with t = (T — T c )/T c . Eqs. © - 
( p3| ) were evaluated previously only for large L. Here 
we present the numerical evaluation of Eqs. ([16]) - p^ ) 
for arbitrary L < 32 without further approximation for 
d = 5 including Wcgner corrections and lattice terms. 
Our strategy of adjusting uq is based on the fact that 
Q at T = T c depends only on uq and that no overall 
adjustment for Q is required since lirnz^oo Q — Qo is 
universal. Thus we adjust uq — 0.93 to the MC data 
H of Q at T c (Fig. 2), then we use the same uq for 
X at T c . For the comparison of x with the MC data 
for Xs at T c we introduce the amplitude A according to 
Xs - AJx/K = A X /(2K C ). Using g X c - 0.1139155 
and adjusting A = 0.678 yields the solid line in Fig. lb. 
At T ^ T c we determine ao = 2.87 from the bulk suscep- 
tibility Xs — 1.322i _1 of series expansion results [ pj| . 

In Figs, lb-3 our analytic result (solid lines) is com- 
pared with the MC data of Ref. J8|. We conclude that 
our one-loop finite-size theory based on H satisfactorily 
describes the existing MC data for 4 < L < 22, both 
at T c and away from T c (Fig. 3). We attribute the re- 
maining deviations of Q for small L to the (expected) 
inaccuracy of our one-loop approximation. At T = T c 
our analytic results approach the lowest-mode results 
liniL^oo XsL- 5/2 = po = 1-757 and Q = 0.4569 (hor- 
izontal lines in Figs, lb and 2) from above, in particular 
our theory predicts a (weak) maximum of XsL~ 5 ^ 2 at T c 
(similar to that in Fig. la for n = oo) that has not yet 
been detected in the MC data ||. Our theory also pre- 
dicts a nonmonotonic L dependence of Q at T c (Fig. 2) 
and of the scaled magnetization < |$| > L 5 / 4 at T c . 

Finally we answer the question to what extent the MC 
data in Figs, lb-3 can be described by the finite-size 
scaling forms of Xscai = ^JXscai and Q sca i derived pre- 
viously (Eqs. (76) - (88) of Ref. @) on the basis of H. 
These scaling forms neglect Wegner corrections and lat- 
tice effects. We have found that the same scaling func- 
tions can be derived on the basis of H provided that a 



3 



smooth cutoff is used. The corresponding scaling func- 
tions depend on the two scaling variables x = t(L/^) 2 
and y = (L/Zn) 4 ~ d where £o oc a is the amplitude 
of the bulk correlation length and Iq oc u^ d ^ is a sec- 
ond reference length. Thus, instead of uq and an, we 
now have Iq and £n as adjustable parameters. Since the 
one-loop results for \ and Xscai differ at 0(u 2 ) one must 
allow for a different amplitude A sca i A in the adjust- 
ment of Xscai to Xs- Using the same strategy of adjust- 
ment as described above we find in = 2.641 from Q at 
T c and A sca i = 1.925 from Xs = A sca iXscai/{2K c ). Fi- 
nally we determine £n = 0.396 from the one-loop bulk 
result Umi-o limz-oo X«< = A scal $/(2K C ) = 1.322. The 
corresponding scaling results are shown in Figs, lb-3 as 
dashed lines. We identify the significant departure of the 
MC data for Xs at T c from the dashed line for L < 12 as 
a lattice effect that is well described by our full one-loop 
theory (solid line in Fig. lb) but which is not captured 
by the scaling form. 



which, together with the fitting formula Eq.(32) of Ref. 
@, implied Co = 0.549 and l Q = 0.603. This formula 
omits the leading Wegner correction oc L( 4_d ) d / 4 and a 
negative lattice term oc L~ d l 2 [compare our Eq.([l3|)] and 
therefore implies an increasing XsL~ 5 ^ 2 (Fig. 9 of Rcf. 
H) towards hin^^oo x s L~ 5 / 2 = po = 1.91, in contrast 
to the decreasing XsL~ 5 ^ 2 with po — 1.76 of our one-loop 
theory. More accurate MC data would be desirable which 
could distinguish between our quantitative predictions in 
Figs, lb and 2 and those implied by the analysis of Rcf. 
||. It would also be desirable to determine £n for the 
d = 5 Ising model (e.g. from series expansion results) in 
order to resolve the disagreement between our prediction 
for £o and that of Ref. ■ 
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FIG. 2. Moment ratio Q at T c for d = 5 and n = 1. 
and dashed lines approach the lowest-mode line for L - 
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FIG. 3. Temperature dependence of susceptibilities for 
d = 5 and n = 1: l(T 2 Xs for L = 4 and l(T 3 Xs for L = 12 . 

This failure of the scaling form for L < 12 was first ob- 
served by Luijten et al. ||. We see, however, that there is 
good agreement of our scaling results with the MC data 
for L > 12, contrary to the disagreement found in Ref. || . 
The latter disagreement is due to the (unjustified) iden- 
tification M J — K,Xs = X corresponding to A sca i = 1 
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